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I. INTRODUCTION 



Military officers, unlike other professionals, cannot practice 
their profession except in time of war. Consequently, throughout 
much of history the military has developed, or has sponsored 
and supported the development of methods and techniques that will 
permit it to practice its profession in time of peace. One such 
technique is based on the computer simulation of war and is 
know as "war gaming”. 

War gaming is a technique which may be used by a force comman- 
der to evaluate the relative merits of various courses of action 
available to him or her for the accomplishment of his or her 
mission. In training situations, war gaming allows the commanders 
and future commanders to exercise and develop their analytical 
skills, to examine the validity of established doctrine, and to 
explore the viability of evolving tactical concepts. The general 
purpose of a war game is an "educational” (decision making experience) 
or an "analytical" (information evaluation experience) type of game. 

The Interactive Simulation of Engagements At Sea (ISEAS) is a 
PC-based naval war game under development at the Naval Postgra- 
duate School. This game permits the evaluation of various tactical 



and 


engineering options 


under 


the 


constraints 


of a specified set 


of 


rules 


( Figure 1.1) . 


As a 


naval 


tactical 


war game, ISEAS must 


provide a 


realistic simulation 


of 


detection , 


engagement, Combat 



Information Center (CIC), and destruction of the threat. 
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Figure 1.1 ISEAS 



The purposes of ISEAS are to give the players experience in 



supervising an air attack and in examining the effects of engineering 
changes to the air defense system. Because it can be played over 
and over again, the war game makes it possible to do what 
cannot be done in the field, that is, to vary characteristics and to 
extend the scope and values of limited peacetime experiments. 

To be of value, ISEAS must provide a realistic radar detection 
model for the airborne threat. This thesis describes the development 
of such a detection model. If the detection of a target cannot 
be satisfactorily modeled, then its engagement and destruction 
cannot be simulated, since in the sequence of events it is nece- 
ssary to detect a target before a weapon system can be brought 
to bear. 
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II. SYSTEM DESCRIPTION 



This report contains a description of the computer model of 
a ship-borne air defense radar system. The model is written in 
C-language for use on a PC microcomputer . A general description 
of the model is presented in Figure 2.1 and is followed by a 
detailed documentation . 



The 


function 


of a modern air 


defense radar system 


is to 


detect , 


track , 


select, designate, 


and direct 


missiles 


and guns 


to an 


intercept 


with unfriendly 


aircraft or 


missiles . 


Such a 


system 


consists 


of five main subsystems: 







- Search radar(s) 

- Tracking radar(s) 

- Data processor(s) 

- Guided missile interceptors 

- Command and control network(s) 

The search radar obtains the initial information about approaching 
air targets. This information is analyzed in the data processor, 
which selects specific targets for tracking and designates the 
tracking radars. The tracking radar then provides the data processor 
with extremely accurate data indicating the target’s location and 
speed. From this information, the data processor chooses an inter- 
cept point for the missile trajectory or projectile to engage 
each target. Guided missiles are then launched by the launching 



equipment and are guided to the target, possibly by the radar 




Figure 2.1 Radar Detection 
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trackers and the data processor. The function of the command and 



control subsystem is to interconnect the various subsystems. 

A. RADAR 

Radar is an electromagnetic system for the detection and loca- 
tion of objects. It operates by transmitting a particular type 
of waveform, such as a single pulse. An elementary form of radar 



consists 


of a transmitting antenna emitting electromagnetic re- 


diation 


generated by an oscillator of some sort, a receiving 


antenna , 


and an energy-detecting device or receiver. A portion 


of the 


transmitted signal is intercepted by a reflecting object 



(the target) and is reradiated in all directions. It is the energy 
reradiated in the direction of the receiving antenna that is of 
prime interest to the radar. The receiving antenna collects the 
returned energy and delivers it to a receiver, where it is pro- 
cessed to detect the presence of the target and to extract it’s 



location 


and possibly it’s relative radial velocity. The distance 


to the 


target for a monostatic radar is determined by measuring 



the time taken for the radar signal to travel to the target and 
back. Radio waves travel at 328 yds every microsecond, i.e every 
millionth of a second, which means that every microsecond represents 



164 yds 


of range. 


The 


echo returned to the radar by a simple, point-source 



target will be an exact reproduction of the transmitted signal 



shifted 


in time by an amount corresponding to the range delay and 


shifted 


in frequency by the "Doppler shift" due to the target 



relative radial velocity. [Ref. 1] 
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If "n" echo pulses, all of the same signal-to-noise ratio, were 
integrated (the process of summing all of the radar echo pulses 
for the purpose of improving detection) by an ideal predetection 
integrator, the resultant signal-to-noise (power) ratio would be 
exactly n times that of a single pulse. If the same n pulses were 
integrated by an ideal postdetection device, the resultant signal- 
to-noise ratio would be less than n times that of a single pulse. 
This loss in integration efficiency is caused by the nonlinear 
action of the second detector. Experimental measurements of the 
operator’s ability to detect signals on an A-scope showed the 
integration-improvement factor to be root (n), or 1.5 dB per dou- 
bling of the number of pulses integrated. 

Since the range to the target can be predicted with reasonable 
confidence, these ranges have been viewed as being determinis- 
tic in character. 

Two basic characteristics of a target can be obtained from its 
radar echo; metric characteristics and signature characteristics . 
Metric characteristics give the target location in terms of range 
and angular measurements. Signature characteristics give the target 
size and shape. 

B. PULSE MODULATED SEARCH RADAR 

A pulse modulated air seach radar is an active system which 
emits spaced pulses of electromagnetic energy. These pulses of 
energy illuminate an airborne target and are reradiated by it. A 
portion of the reflected energy is received by the radar set, 
and the range to the target is determined by the time differential 
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between the transmission of a pulse and the receipt at the radar 



receiver of the target echo. 

The range to the target from the radar is determined by the 
following relation: 

R = C_SPEED * t / 2 (1) 

where 

R = Range to the target 

t = Time difference between the transmission of a 
radar pulse and the receipt of the target echo 
C_SPEED = The Speed of Light 

C. TRACKING RADAR 

A tracking radar is one whose function is to report the 

position of a moving target and to follow its movement in both 
range and angle. The function of these radars is not the detection 
of a target. Their range performance is considered in a different 
sense from that of a surveillance or search radar. Typically, a 
tracking radar acquires new targets on the basis of approximate 
target position information furnished to it by a search radar. 



The range-angle 


gates of 


the 


tracking radar 


are 


positioned at the 


range and angle 


specified 


by 


the 


search 


radar . 


In some cases the 


tracking radar 


is scanned 


in 


both 


range 


and 


angle over a limited 



volume of space surrounding the designated position. If the echo 
from the target appears within the gate with an adequate signal- 
to-noise ratio, the tracking radar "locks on" to the target echo. 
Therefore, the problem of analyzing the range performance of a 
tracking radar is reduced to that of finding the range at which 
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the signal-to-noise ratio has the value required for successful 



target acquisition or lock-on. 

D. TIME-STEP AND EVENT-STORE GAMES 

Computer games are conducted according to either a time-step 
or event-store method. In the time-step method all moves are made 
for a fixed time interval. The computer calculates the positions 
of the threat at regular time steps. At the end of each step, 
it examines each interaction that might have occurred, decides if 
it did occur, and, if it did, determines the outcomes. 

In the event-store or critical event method, the computer, 



instead 


of advancing the 


game 


by constant 


time 


increments, jumps 


ahead to 


the 


next event 


that 


has a probability 


of occurring and 


that is 


vital 


to the outcome of 


the game. 


No positions, distances, 


or times 


are 


determined 


until 


they can be 


used . 





As each new event is generated an "event word" is formed. This 
word tells when the event will occur, what type of action is to 
be taken, and what offensive and defensive units are involved. The 
word is then stored according to its time of occurrence in a 
chronological list of future events. 

START > TAKE NEXT EVENT DETERMINE TYPE EVENT 

/;\ I 

• \ 

f f 

; I 

• V/ 

STORE NEW EVENT OR CARRY OUT ROUTINE FOR 

< 

EVENTS AS REQUIRED THIS TYPE OF EVENT 



Figure 2.2 Event Store Routine 
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The sequence of the game is as follows: Input the data and the 



program into the 


computer. The 


machine 


then 


calculates 


and stores 


in proper 


order 


all the events 


that 


can be 


predicted 


at 


the 


start of 


the game 


The program 


then 


takes 


the first 


event 


in 


store, calculates 


the outcome, 


stores 


this , 


forms words for 


any 


required new events, and removes any 


units 


that will 


no longer 


affect the 


game 


results. With 


further 


events in store, 


they 


are 


considered 


in order, until the 


last event is 


considered 







One advantage of the event-store method is that target detec- 
tions, weapon releases, weapon kills, etc, occur in proper time 
sequence. Another advantage is the ease and flexibility of pro- 
gramming. This is particularly true in the C-language when an 
effort is made in the original coding to (1) allow space in the 
computer for expansion, and (2) to make each event a selfcontained 
routine. Since experience shows that all games need some modifi- 
cation, this feature is well appreciated. One disadvantage of the 
event-store technique is the large amount of memory that must be 
set aside for storing future events. 

E. DETERMINISTIC MODEL 

Computer models, like other game models, may be either determinis- 
tic or stochastic in their treatment of chance events. A deterministic 
model has the property that the result of a modeled event is com- 
pletely determined by the data which describes the situation. In 
such a model, the results of probabilistic events are determined 
by expected values for the events; e.g. if the probability of 
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aircraft survival over a given route is 0.20, then if ten aircraft 



fly along this route, the deterministic model says exactly eight 
are killed and two are alive at the end of the route. A stochastic 
or Monte Carlo model has the property that the result of a modeled 
event is determined by the interaction of the model rules with a 
random process (or processes). For instance, the stochastic model 
might determine how many of the ten aircraft survive by comparing 
ten random numbers with the survival probability 0.20; for each of 
ten numbers randomly drawn between 0.0 and 1.0 less than 0.20 of 
aircraft survives. Repeated runs of a stochastic model reveal the 
distribution of possible outcomes, not just an expected outcome. 
Stochastic models are also called probabilistic models. 
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Ill . 



ENGINEERING DESIGN 



The radar computer program presented here was developed to 
provide the radar performance required in ISEAS. The program 
computes the performance of ship-based search and fire control 
radars against targets in the presence of clutter and jamming. 
The effects of attenuation due to atmosphere and the accuracy 
of radar measurements are included in the program. The program 
uses the basic equations for radar performance described by N. I. 
Skolnik [Ref. 1] and S. A. Hovanessian [Ref. 2], 

The development of a computer program requires an understan- 
ding of the problem, a knowledge of the instructions which the 
computer is capable of performing, the language by which it ope- 
rates, the capacity of the memory, and the method of data 
represen tation . 

The first step in the preparation of the program is the 

contruction of a generalized model, such as the flow chart on 

Figure 3.1. Each of the modules presented on the figure are pre- 

sented below. 

A. RADAR HORIZON 

A radar signal is an electromagnetic wave that travels in 
straight lines. This means that a target can only be detected 
if is on a 1 ine-of-sight from the radar. That is, it will not be 
visible if it is below the horizon. 

The distance from a radar at height h.1 to the horizon (D) may 
be determined by simple geometry. The radius of the earth is 
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Figure 3.1 Flow diagram radar equation 
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Compute "K" 



I 

Compute PSS 



! 

Compute Pc 



f 

Compute 

errors 



_! 

Output 



Figure 3 . 1 



Functions 



Equations Used 



k_f actor () [Eqn . 17] 

n K M = Factor that depend 
of n n H and false 
alarm number 

r^sprob () [Eqn. 16] 

PSS= Single scan 
probability 



r_cprob () 

Pc= Cumulative 
probability 

r_r angler ror () 
r_elv_error () 
r_bear_error () 



[Eqn. 18] 



[Eqn. 21] 
[Eqn. 23] 



pass__value () 

Range 

Elevation 

Bearing 



Flow diagram range equation (cont) 
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approximately 3450 N.M. Because the radar direction changes due 
to refraction of the radar signal by the atmophere, the radar 
earth’s radius is taken as 

4/3 * 3450 = 4600 N.M. 

assuming a factor of 4/3 due to refraction. 

Thus, the distance to the reference radar horizon is given by [Ref. 3] 
D = 1 . 23 root ( h . 1 ) 

where D is in nautical miles and h.1 is in feet. 

If the target is at a height h.2, in feet, the maximum range 

from the radar at which it can be detected (R ) is given 

horiz 

by [Ref. 3] 

R = 1.23 * [ root ( h.1) + root(h.2) ] (2) 

horiz 

For example, if 

h.1 = 30 feet 
h.2 = 400 feet 

R = 1.23 * [ 5.48 + 20.0 ] = 33 N.M. 

horiz 

B. MINIMUM RANGE 

When range measurement is made by using a pulsed waveform 
and by measuring the elapsed time between the transmitted pulse 
and the received signal, there are limitations to the minimum 
range that can be measured. The pulse width is one limitation, 
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since the receiver cannot receive a signal until the transmitter 
is turned off. Another limitation is the time interval required 
to turn the receiver on after the transmitter pulse is turned 
off. The latter is referred to as the recovery time. The minimum 
range is given as [Ref. 4]. 

R = C_SPEED *(T+RT)/2.0 (3) 

min 

Where 

T = Pulse width 
RT = Recovery time 

For example, if 

T =3.0 msec , RT = 1 . 0 msec 

Then 

R = 984 ft / msec * ( 3.0 + 1.0 ) msec / 2.0 
min 

R = 1968 ft. 
min 

A radar system that must have short range capability must use 
a very short pulse width with a quick receiver recovery time. 

C. MAXIMUM UNAMBIGUOUS RANGE 

During one complete transmit- receive cycle, a radar pulse 
must travel from the radar to the target and return. Echo signals 
received after an interval exceeding the pulse repetition period 
are called multiple- time-around echos. They can result in erroneous 
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or confusing range measurements. That is, the radar "sees" un- 
ambiguously only those signals which return to it between successive 
pulses. Consequently, the maximum unambiguous detection range 

(R ) is determined by the pulse repetition frequency ( PRF ) 

unamb 

from the equation [Ref.l] 

R = C_SPEED / (2.0 * PRF) (4) 

unamb 



For example, if 



Then 



PRF = 1000 pulses /sec 



R = 164,000 yds = 82 nautical miles, 

unamb 



One method of distinguishing mult iple-time-around echos from 
unambiguous echoes is to operate with a varying PRF. The echo 
signal from an unambiguous range target will appear at the same 
place on the A-scope on each sweep regardless of whether the PRF 
is modulated or not. 



D. SYSTEM LOSSES 

There are many sources of power loss in the signal and the 
echo. The overall loss is called the system loss and is the sum 
of all of the losses. For that reason, if the effect of the system 
loss is included in the radar range equation, the equation will 
give a more accurate picture of practical achievable ranges. 
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The transmitting losses include such items as power trans- 



mission efficiency, wave guide and antenna losses, and the duplexer 
loss. The total of these losses, usually referred to as plumbing loss, 
in a typical radar may be approximately 1 to 3 decibels. [Ref. 2] 
In the receiving loop, the losses will include the antenna, 
wave guide, RF or parametric amplifier, mixer, and IF amplifier. 
These losses may be on the order of 1 to 3 decibels. 

In addition to the above losses, the electromagnetic signal 
traveling through the atmosphere suffers another loss which is 
termed atmospheric attenuation. This loss is caused primarily by 
absorption of energy by atmospheric gases . For frequencies below 
10 GHz, it is between 0 to 3 decibels depending on the transmitting 
frequency, the antenna elevation angle and the target range. 

As the antenna beam sweeps the target, the antenna gain illu- 
minating the target varies. The radar range equation, however, 
uses the maximum gain of the antenna along the boresight axis as 
the reference. The gain variation of a scanning antenna as it 



illuminates 


the 


target 


will give rise to 


an 


antenna beam 


shape 


loss during 


the 


scan . 


These losses, may 


be 


approximately 


1 to 


3 decibels, 


and 


are 


known as beam shape 


losses . 





E. 


ATMOSPHERIC 


ATTENUATION 










Attenuation 


results when a portion 


of the 


energy incident on 


the 


molecules 


of 


the atmosphere 


is 


absorbed 


as heat and is lost. 


The 


reduction 


in 


radar signal 


power 


when 


propagating over a 

-(& * R/5) 



distance R and back (two-way path) may be expressed as 10 



where & is the (one-way) attenuation coefficient measured in units 



2k 



of distance. At frequencies below about 1 GHz, the effect of 
atmospheric attenuation is negligible. Above 10 GHz, it becomes 
increasingly important. The large attenuations experienced at milli- 
meter wavelengths is one of the reasons long range radars are 
seldom found above 35 GHz. 

From Figure 12.10 of Reference 1, the attenuation coefficient 
for frequencies between 1 GHz and 15 GHz is given as 

& (Db / N.M) = 0.003001 + 0.000393F 

Where F is r_freq / GHz 

F. RECEIVER NOISE 

Noise is the important factor that limits the receiver sen- 
sitivity. It is unwanted electromagnetic energy which interferes 
with the ability of the receiver to detect the echo. No matter 
whether the noise is generated by an internal thermal mechanism 
(noise generated by the thermal motion of the electrons in the 
receiver input stage) or by some other mechanism, the total noise 
at the output of the receiver is the thermal-noise power obtained 
from an ideal receiver multiplied by a factor called the noise 
figure that is a function of a type of receiver. 

Thermal noise = K * T * Bn 
where 

T = Temperature in degrees Kelvin 
K = Boltzmann’s constant 
Bn = Receiver bandwidth in Hz 
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And the total receiver noise ( n_rec ) can be written 



n_rec =F*K*T*Bn (5) 

where F is the Noise figure 



For example, if 

Bn = 5.0 MHz 
F =5.0 



Then 



-23 6 

n_rec = (5.0) * (1.38 * 10 ) * (290) * (5.0 * 10 ) 

-13 

n rec = 1 . 005 * 1 0 W 



G. NUMBER OF PULSES 

The number of integrated pulses per scan (n) used in the equation 
for the probability of detection can be calculated for many regular 
scan patterns. For a fan beam sweeping in azimuth, n is directly 

proportional to the antenna beamwidth area (BA) and PRF, and 
inversely proportional to antenna scan rate. Thus, according to 
Reference 5 

ST * PRF * BA 

n = (6) 

SFOV 

Where 

ST = Scan time = 360 / AZR , sec 
AZR = Antenna azimuth scan rate , deg / sec 
PRF = Pulse repetition frequency , pulses / sec 
BA = Beam area = ( BAZ ) * (BEL), square degrees 
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BAZ = Antenna azimuth beam width, deg 



BEL = Antenna elevation beam width, deg 
SFOV = Search f ield-of-view, square degrees 
= 360 * (BEL) 

For example, if 



Then 



And 



AZR 


= 60 deg/sec 




PRF 


= 240 pulses/sec 




BAZ 


= 2.5 deg. 




BEL 


= 20 deg. 




ST 


= 360 / 60 = 6 sec 




BA 


= (2.5) * (20) = 50 square 


deg . 


SFOV 


= (360) * (20) = 7200 square deg 




(6) * (240) * (50) 




n = 


= 10 


pulses 




7200 




MAXIMUM DETECTION RANGE 




The maximum 


detection range is a 


useful 



determining the maximum range at which an airborne target 

be detected. The maximum detection range (R ) can be giver 

max 

the form [Ref. 1] 



i n 
can 
in 
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( 7 ) 



2 2 

4 Pt * Gr * Lamd * RCS 

R = 

max 3 

( 4*PI ) * Pm 



Where 

R = Maximum detection range, meters 
max 



Pt = Transmitted power, Watts 

Gr = Antenna gain, non-dimensional 

Lamd = Wavelength, meters 

RCS = Target radar cross section, sq . meters 
Pm = Minimum power required at the receiver for detection 
at some desired probability of false alarm. 



For example, if 

Pt = 1 . e6 watts 
Gr = 1000 = 30 db 

Lamd= 1 . 0 meter 
Pm = 5.e-13 watts 
RCS = 1.0 Sq . meter 

Then 

R = 178,000 meters = 90 Nautical miles 
max 

The value of Pm is usually determined in relation to the radar 
receiver noise in the form [Ref. 5] 

Pm = (S/N) * Nr (8) 

c 
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Where 



(S/N) = Signal to noise ratio to give the required 

c probability of false alarm 

Nr = Receiver noise, Watts 



I . ELECTRONIC COUNTER-COUNTERMEASURES 

All radars, civilian and military, must be able to operate in 
the crowded electromagnetic environment that results from the 
transmissions of other radiating sources, both in the radar’s own 
band as well as from other bands. Military radars, however, must 
also operate in a hostile environment where they may be subjected 
to deliberate interference designed to degrade their performance. 
The various methods for interfering electronically with radar are 
called electronic countermeasures , or ECM. Active ECM is sometimes 
referred to as jamming. There is also passive ECM, such as chaff, 
which reflects radar energy to create clutter or a false target. 
The methods employed to combat ECM are called electronic counter 
countermeasures , or ECCM. 

As a general rule, good radar design can reduce vulnerability 
to electronic countermeasures. Good design is based on maximizing 
the ratio of the signal energy to noise power per hertz (E/No), 
as well as employing techniques to reduce mutual interference . 
A radar should be designed more conservatively, with large power 
and large antenna aperture, than the minimum required for marginal 
detection . 
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1 . Noise Jamming 



Internal receiver noise generally limits the sensitivity of 
most microwave radars. Raising the noise level by external means, as 
with a jammer, further degrades the sensitivity of the radar. A 
sensitive receiver with a low noise figure may be desirable in 
many civilian applications of radar, but in military applications 
it is not always an asset since it makes the receiver more 
vulnerable to jamming. 

Noise jamming, which is one of the most effective ECM tech- 
niques, involves generation of noise signals which are spread over 
the frequency bandwidth of the receiving radar. Noise jamming, in 
effect, increases the noise level of the receiving radar, reducing, 
and in many cases eliminating, its capacity to detect and/or 
obtain information on targets. Locating the jammer might not be 
as important as locating the attacking missiles or aircraft 
screened by the jammer, especially when the jamming aircraft is 
beyond the range of defensive weapons. A jammer that operates 
outside the range of normal defenses is known as a stand-off 
jammer. An important ECCM tactic is to engage hostile jammers 
with home-on-jam (HOJ) missile guidance. If the jammer is radia- 
ting, HOJ is generally a better guidance technique than is radar 
guidance . 

The jammer power receiver at the radar varies inversely as the 
square of the distance between radar and jammer (one-way propaga- 
tion), while the radar echo power varies with distance inversely 
as the fourth power (two-way propagation). There will be some 
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distance at which the radar echo will burn through the jammer 



signal. This is also called the self-screening range, or the cross-over 
range . 

2 . Single Scan Probability 

For a jammer that is emitting a power of Pj watts per 
hertz in the frequency band of the radar receiver, the jamming 
signal Nj seen by the radar receiver is [Ref. 5] 

2 

Pj * B r * Gj * Gr * Lamd 

Nj = (9) 

2 2 
( 4*PI ) * Rj 

where 

Pj = Jammer power, watts per hertz 
Br = Radar receiver bandwidth, hertz 
Gj = Jammer gain (antenna focusing) 

Gr = Radar antenna gain 

Lamd = Wavelength, meters 

Rj = Range from jammer to radar, meters 



If the assumption is made that the jammer signal Nj is additive with 
the receiver noise Nr, then the signal to noise ratio in the 
receiver is 
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2 2 

Pt * Gr * Lamd * RCS 

~ 3 4~ 

(4*PI) * Rt 



S/N = 



Nr + Nj 



Pj * Gj * Gr * Lamd * Br 

Nr + 



( 10 ) 



where 



2 2 
(4*PI) * Rj 



Rt = Range from target to radar in meters. 



Defining the reference target range ( Ro ) , and the reference jamming 
range (Rk) as 



2 2 

4 Pt * Gr * Lamd * RCS 



Ro 



(id 



(4*PI) * Nr 



Pj * Gj * Gr * Lamd * Br 



Rk 



( 12 ) 



(4*PI) * Nr 



The characteristics of a given radar system against a target of 
specified radar cross section can be specified by the value of Ro, 
which is the range for a signal to noise ratio of unity, or a S/N of 
zero decibels. 



Then 



(Ro/ Rt) 



S/N = 



1 + (Rk/ Rj) 



(13) 
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J. PROBABILITY OF DETECTION 



The probability of target detection depends on the value of 
the signal-to-noise ratio (with or without jamming) and the threshold 
setting in the detection circuitry. Higher values of S/N signify 
higher values of target return power for a given noise power. 
Higher S/N thresholds for detection raise the level of echo power 
at which target returns will be detected, thereby degrading detection 
capability. But, on the other hand, higher thresholds reduce the rate 
of false alarms (noise crossing the threshold level) making the 
detection of targets more reliable. 

Expressions relating the single scan probability of detection 
(PSS), as a function of the signal-to-noise ratio with the false 
alarm probability, as a parameter are given below. The false alarm 
probability Pn refers to the condition where the noise exceeds 
the threshold setting resulting in a false detection. 

With this information, the designer of a radar system can 
determine the probability of detecting a given target as a func- 
tion of the reference target range (Ro). 

The S/N at any range R can be calculated using the equation 



4 

S/N = (Ro/R) (14) 



and the single scan probability of detection at range R is given 
by [Ref. 5] 



PSS = exp [ -K / (S/N) ] (15) 
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PSS 



(16) 



4 

= exp [ -K * (R/Ro) ] 



where 

R = Range to target 

K = A factor for each radar which is a function 
of n and m 

n = No. of pulses per look (scan) 
m = False alarm number 

The K-factor can be calculated as follows [Ref. 5] 

If 1 . e-6 <= m <= 1 . e-1 2 

K = 10 log K 
db 

K = -6 log n + .4 log m + 7.2 
db 

■ 1 ( K ) 

K = 10 db (17) 

For example, if 

For n = 10 , m = 10e-12 , R = Ro 

K = -6 + 4.8 + 7.2 = 6 
db 

K =4. 

Thus, the single scan probability of detection becomes 

4 

PSS = exp [ -4. ( 1 ) ] = .018 
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and for (R / Ro) = 1/2 



PSS = .78 



For n = 10 , R = Ro and various false alarm numbers 



log m = 2 4 6 8 10 12 

PSS = .36 .22 .105 .068 .036 .018 

Given the probability of detection in a single scan of the 

target, the cumulative probability of detection over several scans 

can be calculated. In scanning radars, targets are illuminated 

every scan or radar frame. For the "s" th scan of the target, there 

s 

will be a certain probability of detection ( P ) based on the target 

d 

range and radar cross section, both of which can vary with each 

_s 

scan. Thus, P is generally a function of time. The cumulative 

d 

probability of detection is obtained by "adding” the probability 

of detection at each scan. The term "adding” is used in a pro- 
bability sense rather than an algebraic sense. The probability 

that the target has been detected after "s" scans, is given by 



_ s 1 2 s 

P = 1 - (1 - PSS) * (1 - PSS) ...(1 - PSS) 

d 



If PSS is essentially a constant value 
can be simplified as [Ref. 6] 



s 

for each scan, then, P 

d 



P 



s 



d 



s 

(1 - PSS) 



(18) 
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R refers to the range at which a 90 % cumulative probability of 
90 

detection is achieved. 

The cumulative probability has sometimes been proposed as a 
measure of the detectability of a radar rather than the single 
scan probability of detection. In practice, the use of the cumu- 
lative probability is not easy to apply. Radar operators do not 
usually use such a criterion for reporting detections. They seldom 
report a detection the first time it is observed, as is implied 
in the definition of cumulative probability. Instead, the criterion 
for reporting a detection might be a threshold crossing on two 
successive scans, or a threshold crossing two out of three scans, 
etc. In track-while scan radars, the measure of performance might 
be the probability of initiating a target track rather than a 
criterion based on detection alone. 

K. CW RADAR 

A CW radar to measure both range and velocity of a single 
target is accomplished by broadening the transmitted frequecy 
spectrum by frequency modulating the carrier frequency. Range, in 
terms of time, can be determined by measuring the difference between 
the received and transmitted frequencies. 

CW radars are used also for CW missile guidance systems and 
for measuring velocity. For the latter, The transmitter generates a 
continuous (unmodulated) oscillation of frequency ( f o ) , which is 
radiated by the antenna. If the target is in motion with a velocity 
Vr relative to the radar, the received signal will be shifted in 
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frequency from the transmitted frequency fo by an amount +- fd 



given by [Ref. 1] 



2 * Vr 

fd = 

Lamd 

The received echo signal at a frequency fo +- fd enters the radar 
via the antenna and is heterodyned in the detector (mixer) with a 
portion of the transmitter signal fo to produce a doppler frequency 
shift, fd, that gives the relative velocity of the target. 

CW radars are smaller than corresponding pulse radars, although 
two antennas are required, and can operate against a target to 
almost zero range. As in surveillance radars, tracking radars 



are 


also 


limited 


by 


coverage 


volumes and 


target 


cross sections. 


The 


same 


considerations apply 


as in the 


case of 


search radars. 


The 


target 


must 


be 


within the 


i coverage 


limits 


and within those 



limits 



0 . 25 

R = Ro [ RCS ] 



(19) 



where 



Ro = R for RCS equal to 1.0 square meter 

But for tracking radars we will operate with relatively narrow 
beamwidths at mid-frequencies (2000 to 6000 MHz). 
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L. ACCURACY OF RADAR MEASUREMENTS 



The parameters of interest in radar applications are the range 
(time delay) and the angle of arrival. Therefore, the theoretical 
accuracies of radar measurements are given below. 

1 . Range Accuracy 

The measurement of range is the measurement of time delay 



Tr = 2 * R / C_SPEED 

In terms of the effective bandwidth, the time-delay error (&Tr) 
becomes [Ref. 1] 



1 

&T r = (20) 

0.5 

betha * (2 * S/N ) 



where betha, the effective bandwidth, is 
2 * Bn 

betha = 



T 

and the range error ( &R ) will be given as 



&R = (C_SPEED / 2) * &T r (21) 

The rms time-delay error of a waveform with a uniform frequency 
spectrum of width B is given as [Ref. 1] 

1 . 752 sin x 

&Tr = * (22) 

0.5 

PI * B * (2*S/N) x 



where the waveform is of the form sin(x) / x , and x = PI * B /freq. 
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2 . Angular Accuracy 



The 


measurement 


of 


angular position is the 


measurement of 


angle 


of arrival 


of 


the equiphase wavefront 


of the echo 



signal. The theoretical error of the angle measurement (&tetha) 
may be derived in a manner similar to the derivations of time error 
(range). For simplicity, the angular error in one coordinate plane 
only will be considered. The analysis will be extended to include 
angular errors in both planes. In terms of the half-power beamwidth 



0.628 * tetha 

&tetha = (23) 

0.5 

(2 * S/N) 

where tetha is the half-power beamwidth in radians, and the units of 
&tetha are the same as those of tetha. [Ref.1] 

3 . Radar Errors Distribution 

When a radar system tracks an aircraft, the target coor- 
dinates measured by the radar, such as range and angle will 
typically fluctuate about the true value. This error can be 

represented by a normal distribution N(m,a), where m is the mean 
2 

or bias, and a is the variance. [ Ref. 6] 

2 

Thus, given X = N(0,1), we can obtain X’ = N(m,a ) by setting 



X 1 = m + (a * X) (24) 
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The C-language library can be used to generate the random number 



and the polar method in Reference 7, to generated standard random 
variables N( 0 , 1 ) 

1. Generate U1 and U2 as U(0,l) , using the function drand48( ) , 

let VI = 2U1 - 1.0 and V2 = 2U2 - 1.0, then 

2 2 
W = VI + V2 

2. If W > 1.0, go back to sted 1. Otherwise, let 

0.5 

Y = [ (-2 Ln W) / W ] 

XI = VI * Y, and X2 = V2 * Y. Then XI and X2 are N(0,1) 
random variables. [Ref. 7] 
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IV. SUMMARY 



From the standpoint of a clear exposition , it is rather unfortunate 
that a true understanding of the radar detection problem is wrapped 
in complexities of statistical theory which do not convey to the 
practicing designer a real feel for the problem. I have attempted 
to alleviate this problem by confining some of the more detailed 
mathematical derivations to simple assumptions that limit the 
validity of the results. 

The chief advantage of ISEAS is replication. Within a reasonable 
amount of time, a game can be replicated sufficiently often to obtain 
a detection which gives an indication of how bad the situation 
might be, how good it might be, and what action could be taken. 
Then, one or more parameters can be varied, another series 
played, and new results obtained. A comparison of the ranges 
reveals the effects of the parameter variation or variations. 

In conclusion, the effectiveness of the model depends on its 
ability to be sensitive to changes. A very sophisticated model 
may take into account many elements of the encounter but in so 
doing may require a large amount of input data, and excessive 
computer run times. In addition, the increased amount of detail 
in the radar model may not yield a significant increase in sen- 
sitivity or increase the validity of the complete simulation. 
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APPENDIX A 



THE RADAR OF ISEAS 









LTjg CLAUDIO MONTERO 









^include "stdio . h M 
^include "math.h" 

^include <rad_head.h> 

float TR, RANGE, ELV, BEAR ; 

int DET; 

y************-****************************************************--*** 
External variables to pass to the calling routines 
RANGE = Range from the radar to the target, 

ELV = Angle of elevation of the radar, 

BEAR = Bearing of the radar (relative to the ship heading) 

**********************************************■*•*■*-*■*-*■*■**•*■*■*■*-*■*-*■*■**•■*•* ^ 

/* Default values for the radar */ 

float PRF = 1 000 . 0 , AA=30 , T=2e-6 , RT = 0 . 5e-6 , EA = 0 . 8, FF = 9.0e9, 

PL=2 . 0 , BL=3 . 0 , BW=5 . 0e6 , ST=290.0, NF=5.0, 

P0WER=3 . 0e5 , AZR=120.0, BAZ=1.5, BEL=10.0, FAN=1.0e10, 

NS=3 , PJ=0.1, GJ J= 0.0, JR=45 . 0 ; 



/* Default values for the target */ 
float RCS=6 . 0 , x_t=70000, y_t=20000, z_t=10000; 
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INPUT DATA DEFINITION: 

PRF = Pulse repetition frequency, AA = Antenna altitude 
T = Pulse width, RT = Recovery time, EA = Effective area 
FF = Frequency, PL = Plumbing loss, BL = Beam loss 
BW = Receiver bandwidth, ST = System temperature 
NF = Noise figure, POWER = Transmitter power, 

RCS = Target RCS, AZR= Antenna scan rate, 

BAZ= Azimuth beam width, BEL= Elevation beam width 
FAN= False alarm number, NS = Number of scans, 

PJ = Jammer power, GJJ = Jammer antenna gain, 

JR= Range from jammer to radar, x_t= X-position of target, 
y_t= Y-position of target, z_t= Z-position of target (meters) 

*****•*•*■*•****•*■*•***■*******•*******■****#*******-****•**********■********** j 



main( ) 

{ 

float r_un_range( ) , radar_hor iz( ) , r_min_rang( ) , r_cprob(), 
t_range( ) ; 

t_range( x_t , y_t , z_t ) ; 

DET = r_un_range( PRF ) ; 

DET += radar_horiz( AA , z_t ) ; 

DET + = r_min_range( T , RT ) ; 

DET += r_cprob(NS); 
pass_value( x_t , y_t, z_t); 
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y****.****.****.*********-**-*******-****#-*****-*********-**************-#-*-* 



Use to calculate the true range to the target in nautical miles 
given the true target position in polar coordinates 

w********************************************************-********-**/ 

float t_range (x_target, y_target, z_target) 
float x_target, y_target, z_target; 

{ 

float TR1 , TR2; 

TR1 = pow(x_target , 2 . 0 ) + pow( y_target , 2 . 0 ) + pow( z_target , 2 . 0 ) ; 
TR2 = pow( TR1 ,0.5); 

TR = TR2 / 1609.0; /* From meters to N . M */ 

return ( TR ) ; 

} 

^w************************************************-********-*-********* 

Use to calculate the radar unambiguous range of the target 
and return a value depending if the target is into the 
unambiguous range or not. 

■a***************-*-**-***-**-*-***-***********-******-*-********-*-****-***-***-*-**/ 

float r_un_range( r_prf ) 
float r_prf ; 

{ 

float RU1 , RU ; 

RU 1 = C_SPEED/(2.0 * r_prf);/* C_$PEED is defined in rad_head.h */ 
RU = RU 1 / 1609.0 ; /* Meter to Nautical Miles */ 

if ( RU >= TR ) 

return (DETECTED) ; 
else 

return (UNDETECTED) ; 

} 



kk 



/— 



Function used to calculate the horizon of radar in N.M 
depending on the altitude of the antenna and the target, 
then return a value depending if the altitude of the target 
allows the radar to detect it or not. 






float radar_horiz ( r_ant, t_alt ) 
float r_ant, t_alt; 



float 


RH, theta, alpha: 






RH = 1 


.23 * ( sqrt( r_ant ) ) ; 


/* 


RH is in N.M */ 


theta 


= r_ant / (RH * 5280.0); 


/* 


Mile to feet * 5280 */ 


alpha 


= t_alt / (TR * 5280.0); 










tetha and alpha will be proportional to the angle of sight 
of the antenna or target respectively. 



if (alpha > theta ) 

return ( DETECTED ); 
else 

return ( UNDETECTED ); 
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/********************************-*********-**-**-**-*-****-***-*-***-******* 
Use to calculate the minimun range for detection, also return 

a value depending if the range is in or out of the minimum 
range zone. 

float r_min_rang ( r_tau , r_taur ) 

float r_tau , r_taur ; /* r_tau = Pulse width */ 

{ /* r_taur = Receiver recovery time */ 

float MR1 .MR; 

MR1 = ( C__SPEED * (r_tau + r_taur)) / 2.0 ; 

MR = MR1 / 1609.0; /* Meters to Nautical Miles */ 

if ( TR >= MR ) 

return (DETECTED) ; 

else 

return (UNDETECTED); 

} 

^**********-*******-************************** 

Antenna gain ( r_gain ) should be used i 
since it includes the losses introduced 
Assuming antenna efficiency is unity 
Ae = effective area 
A = actual antenna area 
regain = 4 * PI * Ae / SQUARE ( Lamd ) 
*************-**************************************************.*.**.*/ 



*#■##■*•*•#■•*■*•*■*•*•*■*-*•*-*■*•*■*•*-*■* 
n the radar equation 
by the antenna. 

Ae = A. 
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float regain ( r_area, r_freq ) 
float r_area, r_freq; 

{ 

float lamda, AG ; 

lamda = C_SPEED / r_freq ; /* given in meters */ 

AG = (4.0 * PI * r_area ) / pow( lamda , 2 . 0 ) ; 
return ( AG ) ; 

) 

The system loss is the sum of all the losses. They will be 
given in Db, but the losses should be in decimals to be used 
in the range equation. 

Plumbing loss r_plumb 

Beam shape loss r_beam 

**********************************************■*■*-*■*-*•*■*-*•*•#■*•»-*■*■*■**■■***■*/ 



float r_losses ( r_plumb, r_beam ) 
float r_plumb, r_beam; 

( 

float RL1, RL2, RL; /* RL1 will be the losses in Db */ 

RL 1 = r_plumb + r_beam; 

RL2 = RL1 / 10 ; 

RL = pow( 1 0 . 0 , RL2 ) ; /* Range equation works with decimal */ 

return ( RL ) ; 

} 
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***W****WW*-****************-***-**W*********-W****-***************--**-** 



The reduction in radar signal power when propagating over a 
distance R and back (two-way path) may be expressed as 
-(& * R/5) 

10 , where & is the attenuation coefficient 

*************-***************#-***#*************-****-***-#-***w-w**#-***ww j 

float r_attenuat( r_f req ) 
float r_freq; 

{ 

float FI, ATT1 , ATT2 , ATT; 

FI = r_f req / 1 . 0e9 ; 

ATT 1 = (0.000393 * FI) + 0.003001; 

ATT2 = -( ATT1 * TR / 5.0); 

ATT = pOw( 10. 0, ATT2) ; 
return( ATT ) ; 

) 

float j_attenuat( j_range , r_f req ) 
float j_range, r_freq; 

{ 

float FI, AJT1 , AJT2 , AJT ; 

FI = r_f req / 1 . 0e9 ; 

AJT1 = (0.000393 *F1 ) + 0.003001; 

AJT2 = -( AJT1 * j_range / 10.0); 

AJT = pOw( 10.0, A JT2 ) ; 
return (AJT); 

} 
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The receiver noise generated by the thermal motion of the 
conduction electrons at the receiver, are given by 
n_rec = K * TS * Br * Fn 

where K is Boltzman’s constant ( 1 . 38e-23 j oules/ degree ) 
TS is System temperature ( 290 K ) 

Br is Receiver noise bandwidth (hertz) 

Fn is Receiver noise figure ( equal to 5 ) 



****************************■#■**■#■*■#■#■*■*■#■#■*■#■*■#*■■*•#■#■*■#■*•*•#■#■#•#■#•*•#■#■*•#■#■*■**•■**■-* 



/ 



float n_rec ( r_bandwid, r_temp , r_figure ) 
float r_bandwid, r_temp , r_figure ; 

{ 

float RN; 

RN =( K_BOLTZ * r_temp * r_bandwid * r_figure); 
return (RN ) ; 

} 

^w*************************-***********-****-*******-******-**-*********- 

Use to calculate the number of pulses per scan (n), that 
will be used in the equation to calculate the single scan 
probability. ( see pag 26 for reference ) 
*************************************•***************•******■*■*****-***•/ 
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float r_pu lses ( r_prf, r_azr, r_baz, r_bel ) 
float r_prf, r_azr, r_baz, r_bel; 

{ 

float r_ st, r_ba, r_sfov, NP; 
r_ st = 360.0 / r_azr; 
r_ba = r_baz * r_bel ; 
r_sfov = 360.0 * r__bel ; 

NP = ( r_ st * r_prf * r_ba ) / r_sfov ; 
return (NP); 

) 

/**********************************■#•***•*•#•#-#■*•*■*■*•*•**■**■*•*•*■*■*•*•*•*■*■*■*■*■*■** 

Radar-target range (Ro) is the characterist ic of a given 

radar system against a target of specified radar cross section 

**********************#*********•*■**■***■**■•*■**■******•#■*■*•*■*■*■*-*•*■*■*•*•*-*■**•*■*• y 

float r_tar_range (r_power, r_freq, t_rcs) 

float r_power, r_f req , t_rcs; 

{ 

float RT 1 , RT2, lamda, AG, RL, RN, ATT, RTR ; 
float r_gain(), r_losses(), n_rec(), r_attenuat( ) ; 
lamda = C_SPEED / r_f req ; 

AG = i^_gain ( EA , FF ) ; 

RL = r_losses( PL , BL ) ; 

RN = n_rec(BW, ST, NF) ; 

ATT = r_attenuat ( FF ) ; 

RT1 = (r_power * ( pow( AG , 2 . 0 ) ) * ( pow( lamda , 2 . 0 ) ) * t_rcs * ATT); 
RT2 = pow( ( 4*PI ) , 3 . 0 ) * RL * RN ; 

RTR = ( pow( (RT1/RT2) , 0. 25) ) / 1 609 . 0; 
return (RTR); 

} 
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y********************************************-*-*-*-*-*-*-*-*-*-*-*-*-****-*-*-#-*-** 



The reference jamming range that will be use to get the 
single scan probability in the presence of jamming. 

float j am_range( j_power , g_jammer, r_f req , r_bandwid) 
float j_power, g_jammer, r_f req , r_bandwid; 

{ 

float lamda , GJ , RJ1 , RJ2 , AG, RL, RN, A JT , RJ ; 
float r_gain(), r_losses(), n_rec(), j_attenuat( ) ; 

AG = r_gain ( EA , FF ) ; 

RL = r_losses(PL, BL); 

RN = n_rec( BW , ST , NF ) ; 

A JT = j_attenuat( JR , FF ) ; 
lamda = C_SPEED / r_freq; 

GJ = pow( 1 0 . 0 , (g_jammer/ 1 0 . 0 ) ) ; /* change Db to decimal */ 

RJ1 = ( ( j_power/1 . 0e9) * AG * GJ * pow( lamda , 2 . 0 ) * r_bandwid * AJT ) ; 
RJ2 = pow((4*PI),2.0) * RN * (RL/2.0) ; 

RJ = pow( (RJ1/RJ2), 0.5)/1609.0; 
return ( R J ) ; 

} 



51 



y*****************************************-*-*-*-**-**--**-***--**-*-****--****-**- 

The signal to noise ratio in the receiver 

************W**W**W***W*W***W****-******WWW*W**WWWWW*W*WW*WW***-****»/ 



float r_sig_noise( j_range) 
float j_range; 

{ 

float RSN1 , RSN2, RTR , R J , RSN; 
float r_tar_range( ) , jam_range(); 

RTR = r_tar_range( POWER, FF, RCS) ; 

RJ = jam_range( PJ , GJJ , FF , BW) ; 

RSN1 = pow((RTR/TR),4.0); 

RSN2 = 1 +( pow( ( RJ/ j_range ) , 2 . 0 ) ) ; 

RSN = (RSN1/RSN2); 
return( RSN) ; 

> 

y****WW#-*W**W-#****-*W***-#W**###*-*WWWWW*-*--**W#*W**#*WW****W-*******-**W* 

K is a factor of the radar equation which is a function of 

n (pulses per scan) and m (false alarm number), where 

K = 10 log K 
db 

K = -6 log n + .4 log + 7.2 
db 

•1(K ) 

K_factor =10 db 
( see example page 34 for reference ) 



*•***•**■**■*■#•**■ 
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float k_factor ( r_false ) 
float r_false ; 

{ 

float k_db1 , NP, KF ; 
float r_pulses( ) ; 

NP = r_pulses( PRF, AZR, BAZ, BEL) ; 

k_db1 = ((-6.0 * log10(NP)) + (0.4 * log 1 0( r_f alse ) ) 

+ 7.2 ) * 0. 1 ; 

KF = pow( 1 0 . 0 , k_db 1 ) ; 
return (KF); 

> 

Single scan probability of detection at range R, is given as 
PSS = exp [ -k /(S/N)] 

( see example page 34 for reference ) 



**■**■**-*■*■■#**■#***#■#**■#*■#*****#**■************************-*******-**-***■*- y 

float r_sprob ( ) 

{ 

float PSS1 , KF, RSN , PSS; 
float k_factor(), r_sig_noise( ) ; 

KF = k_factor( FAN) ; 

RSN = r_sig_noise( JR ) ; 

PSS1 = -( KF / RSN ); 

PSS = exp ( PSS1 ) ; 
return (PSS); 

> 
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/***##***********•**■**•***•**•*■*■**•**■**■*•*■***■**-****••****■*■*•**■•*■*■**■*■*■*■**■**-■**■*■ 
The probability that the target has been detected after "s" 

scans, if we consider PSS constant, is given by 
_s s 

P = 1 - (1 - PSS ) 
d 

( see page 35 for reference ) 

then returns a value depending if the cumulative probability 
is larger than 90 % of detected probability or not. 

**************-******-#**-#******-**W-#****************-**-#*-#*-#**-*-#*-**--H--#-#*y 

float r_cprob ( r_scan ) 
float r_scan; 

{ 



float 


PCI, PC2 , 


float 


r_sprob( ) 


PSS = 


r_sprob( ) 


PCI = 


(1 - PSS 


PC2 = 


pow( PCI , r 


PC = 


(1 - PC2 ) 


if (PC >= 0.90) 



return (DETECTED); 
else 

return (UNDETECTED); 
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This function will calculate range, elevation and bearing 
if the target is detected by the radar, including the range 
and angular error. 

***********************•##■■*■***■##■*■******■*•#■*#■*•*■■#*■#■**•*#■#■****•*#•*•****■*■**■■* y 

pass_value( x_target , y_target, z_target) 

float x_target, y_target, z_target; 

( 

float x_r, y_r , z j, rangel , range2, ran_error, 
elv_error, bear_error; 

float r_rang_error ( ) , r_elv_error( ) , r_bear_error ( ) ; 
if (DET == 4.0) 

{ 

x_r = x_target; 
y_r = y_target; 
z_r = z_target; 

rangel = pow(x_r,2.0) + pow(y_r,2.0) + pow( z_r , 2 . 0 ) ; 

range2 = pow( rangel , 0 . 5 ) ; 

ran_error = r_rang_error ( BW , T ) ; 

elv_error = r_elv_error ( BEL ) ; 

bear^error = r_bear_error ( BAZ ) ; 

RANGE = (range2 + ran_error) / 1609.0 ; 

ELV = ( (asin(z_r/range2) ) * 180.0/PI) + elv_error ; 

BEAR = ( ( atan2(x_r , y_r ) ) * 180.0/PI) + bear_error ; 

} 

printf ( " \nRange is = %.2f N . M M , RANGE ) ; 
printf ( " \nElevation is = $.2f Deg", ELV); 
printf ( " \nBearing is = %.2f Deg " , BEAR ) ; 
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This function will calculate the range and angular error 
as in the formula on page 38. 

**************************************************************************** 



float r_rang_error (r_bandwid, r_tau ) 
float r_bandwid, r_tau; 

{ 



float RE1, RE2 , RER, RSN , RAX; 
float random( ) , r_sig_noise( ) ; 
if (DET == 4) 

{ 

RAX = random( ) ; 

RSN = r_sig_noise( JR ) ; 

RE1 = r_tau / (4.0 * r_bandwid * RSN); 
RE2 = pow( RE 1 ,0.5); 

RER = RAX * RE2 * C_SPEED / 2.0; 
return( RER ) ; 

} 
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float r_elv_error ( r_bel ) 
float r_bel; 

{ 

float RAX, ELER, RSN ; 

float random(), r_sig_noise( ) ; 

if (DET == 4) 

{ 

RAX = random( ) ; 

RSN = r_sig_noise( JR) ; 

ELER = RAX * 0.628 * r_bel / (2.0 * pow((2.0 
return ( ELER ) ; 

) 

} 

float r_bear_error (r_baz) 
float r_baz; 

{ 

float RAX, AZER, RSN; 

float random(), r_sig_noise( ) ; 

if (DET == 4) 

( 

RAX = random(); 

RSN = r_sig_noise( JR) ; 

AZER = RAX * 0.628 * r_baz / (2.0 * pow((2.0 
return ( AZER) ; 

) 

) 



* RSN), 0.5)) 



* RSN), 0.5)) 
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Creates a random variable number N(0, 1 ), using the C-language 
library . 

****#**-*-****-***W-**-K-**##*-***-*«^-**#*-*-*******M-»##-**#**-******W#***W***'-* 



float random( ) 

{ 



} 



double x, 


C\J 

> 

> 

N 

> 


*St 

C\J 

n 

5 


float 


RAX 






while 


( w 


>1.0) 




{ 










x = 


drand48( ) ; 






y = 


drand48( ) ; 






V 1 


- (2 * x) - 1 


• 0; 




v2 


= (2 * y) - 1 


. 0; 




w = 


pow( vl ,2.0) 


+ pow( v2 , 2 . 0 ) 


> 








Z = ( 


(-2. 


0 ) * log( w) ) 


/ w; 


RAX = 


( pow( 2,0.5)) * V 


1 ; 



return( RAX ) ; 
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APPENDIX B 



TRACKING RADAR OF ISEAS 
LTJG CLAUDIO MONTERO 



^include 


"stdio. h" 


^include 


"math . h" 


^include 


<rad_head . h> 


float TR, 


TR_RANGE , T_ELV , T_BEAR ; 


int TRACK; 



EXTERNAL VARIABLES 

TR_RANGE = Range to the target, T_ELV = Angle of elevation, 
T_BEAR = Bearing relative to the ship head. 

/* Default values for the tracking radar */ 
float AAT = 40.0, EAT = 0.6, BWT=4.0e4, CST=4.0, P0WT=1.2e5, 

FFT = 2 . 0e9 , GLT=1.5, BLT = 0.5, ST = 290, NF=5.0, TAZ = 1.0, 

TBEL= 1 .5; 

/* Default values for the target */ 

float RCS= 1 . 0 , x_t=30000, y_t= 20000 , z_t=400; 

/******************************************-*-*•*■*•*■*•*•*•*■*•*’*■-**■*•*•■*•*•*•*-***■-* 
INPUT DATA DEFINITION 

AAT = Antenna altitude (feet), EAT = Effective area antenna 
BWT = Receiver bandwidth, CST = Critical S/N ratio, 

POWT = Transmitter power, FFT = Frequency tracking radar 
GLT = Track gate losses, BLT = Beam losses 

TAZ = Azimuth beam width ( tracker ) , TBEL = Elevation beam width 
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main( ) 



{ 

float t_range(), track_horiz( ) , track_horiz( ) ; 
t_range( x_t , y__t , z_t ) ; 

TRACK = track_hor iz( A AT , z_t ) ; 

TRACK += check_tracking( ) ; 
pass_track ( x_t , y_t , z_t ) ; 

} 

/*#*****************************************-*•*■*■*-*■*■*■**■■**■-*-*-*-*-*-**#■-**■** 
Use to calculate the range to the target in N.M 



float t_range (x_target, y_target, z_target) 
float x^target, y_target, z_target; 

{ 

float TR1 , TR2 ; 

TR1 = pow(x_target, 2. 0) + pow( y_target , 2 . 0 ) + pow( z_target , 2 . 0 ) ; 
TR2 = pow( TR1 ,0.5); 

TR = TR2 / 1609.0; 
return(TR) ; 

) 
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Function used to calculate the horizon of the radar in nautical 
miles depending on the altitude of the antenna and altitude of 
the target. Then returns a value depending on if the radar 
can detect the target or not. 

float track_horiz( truant, z_target ) 
float tr_ant, z_target; 

{ 

float TH, theta, alpha; 

TH = 1.23 * ( sqrt ( tr__ant ) ) ; 
theta = tr_ant / (TH * 5280.0); 
alpha = z_target / (TR * 5280.0); 



if (alpha 


> 


theta ) 


return 


( 


DETECTED ) 


else 






return 


( 


UNDETECTED 



) 

Antenna gain will be used in the radar maximum range equation 

float tr_gain( tr^area, tr_freq) 
float tr_area , tr_f req ; 

{ 

float lamda, AGT ; 

lamda = C_SPEED / tr_freq ; 

AGT = (4.0 * PI * tr_area) / pow (lamda, 2.0); 
return ( AGT ) ; 

} 
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System losses 

Track gate losses tr_gate 

Beam shape loss tr_beam 

■h..*.*.*.*.*.*.*.*.*.*.*.*.*.*.**.***********##***************************-*-*-*-*-*-*-**-*--*/ 

float tr_losses( tr_gate, tr_beam) 
float tr_gate, tr_beam; 

{ 

float TRL1 , TRL ; 

TRL1 = (tr_gate + tr_beam) / 10.0; 

TRL = pow( 10.0, TRL 1 ); 
returnf TRL ) ; 

) 

y*****************************************-*#--*-*-*-*-**-*-*-*-*-*-*-*-*-*-*-*-*-**-***- 
Atmospheric attenuation 

float tr_attenuat( tr_f req ) 
float tr_freq; 

{ 

float FI, ATTR1, ATTR2 , ATTR ; 

FI = tr_freq / 1.0e9; 

ATTR1 = (0.000393 * FI) + 0.003001; 

ATTR2 = -( ATTR1 * TR / 5.0); 

ATTR = pow( 10.0, ATTR2 ) ; 
return( ATTR) ; 

} 
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y******************************************.#.*.*.*.*.**.*.*.*.**..H..#.*.*.*.*.*.*.#.«..**. 



The receiver noise 

float noise_rec( tr_bandwid, r_temp , r_figure) 
float tr_bandwid, r_temp , r_figure; 

{ 

float TRN ; 

TRN = (K_BOLTZ * r_temp * tr_bandwid * r_figure); 
return(TRN) ; 

} 

^***************************************#***-*-*-*-*-*-*-*-*-*-*-**--*-**-*-*-*-H--*-*-** 

Target reference range (Ro) , taking the radar cross section 
equal to one. 

************************************************•*•*•*■**•■*•*•*■**■■*■*•*■*•**•*■*■*/ 
float tr_ref_range( tr_power , tr_f req ) 

float tr_power, tr_freq; 

{ 

float lamda , AGT , TRL , TRN, ATTR , TRR1 , TRR2 , TRR ; 
float tr_gain(), tr_losses(), noise_rec(), tr_attenuat ( ) ; 
lamda = C_SPEED / tr_freq; 

AGT = tr_gain ( EAT , FFT ) ; 

TRL = tr_losses( GLT , BLT ) ; 

TRN = noise_rec(BWT, ST, NF) ; 

ATTR = tr__attenuat( FFT ) ; 

TRR 1 = (tr_power * pow( AGT , 2 . 0 ) ) * pow( lamda , 2 . 0 ) * 1.0 * ATTR; 
TRR2 = pow((4*Pl), 3.0) * TRL * TRN; 

TRR = ( pow( ( TRR 1 / TRR2),0.25)) / 1609.0; 
return( TRR ) ; 

) 
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■* ■* * *■ ■* * * 1 



The maximum range 






/ 



float tr_max_range( t_rcs ) 
float t_rcs; 

{ 

float TRM , TRR ; 
float tr_ref_range( ) ; 

TRR = tr__ref_range( POWT , FFT ) ; 
TRM = TRR * pow(t_rcs, 0.25); 
return(TRM) ; 

) 



check_tracking( ) 

{ 

float TRM; 

float tr_max_range( ) ; 

TRM = tr_max_range(RCS) ; 
if ((0.9 * TRM) >= TR) 
return ( DETECTED ) ; 
else 

return ( UNDETECTED ) ; 

} 
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Calculate the signal to noise ratio 

float tr_sig_noise( ) 

{ 

float TRM, TSN ; 
float tr_max_range( ) ; 

TRM = tr_max_range( RCS ) ; 

TSN = pow( ( TRM / TR ) , 4 . 0 ) ; 
return ( TSN ) ; 

> 

y*******************************************-**--*-*-*-*-**--**-*-*-***--**-*.*--**--** 
Create a random variable number N(0,1), using the C-language 
library . 

****************************************■**■***■*■*■**•■*■*■**■■*•*■*•**••**■■*■*■**■*■■**• ^ 

float random( ) 

{ 

double x, y, z, vl, v2 , w=2.0; 

float RAN; 

while ( w > 1.0 ) 

{ 

x = drand48( ) ; 
y = drand48( ) ; 
vl = (2 * x) - 1 . 0; 
v2 = (2 * y) - 1.0; 
w = pow(v1,2.0) + pow(v2,2.0); 

) 

z = ((-2.0) * log(w)) / w; 

RAN = (pow(z,0.5)) * vl ; 
return( RAN ) ; 
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This function will calculate range, elevation and bearing; 

if the target is detected by the radar 

*************#************•**■•*•**•*******•**•**•*•■**•**•■**■*••***■****■**•**•*•**••*•** ^ 
pass_track(x_target , y_target , z_target ) 

float x_target, y_target, z_target; 

{ 

float x_r, y_r, z_r , rangel , range2, ran_error, elv_error, 
bear_error ; 

float tr_rang_error ( ) , tr_el v_error( ) , tr_bear_error( ) ; 
if (TRACK == 2.0) 

( 

x_r = x_target; 
y_r = y_target; 
z_r = z_target; 

rangel = pow(x_r,2.0) + pow(y_r,2.0) + pow( z_r , 2 . 0 ) ; 

range2 = pow( rangel , 0 . 5 ) ; 

ran_error = tr_rang_error ( BWT , FFT ) ; 

elv_error = tr_elv_error( TBEL ) ; 

bear_error = tr_bear_error( TAZ ) ; 

TR_RANGE = ( range2 + ran_error) / 1609.0; 

T_ELV = ( ( asin ( z_r/range2 ) ) * 180.0/PI) + elv_error ; 

T_BEAR = ( ( atan2(x_r , y_r ) ) * 180.0/PI) + bear_error ; 

) 

printf( " \n Range is = %.2f N . M" , TR_RANGE ) ; 
printf ( "\nElevation is = %.2f Deg " , T_ELV ) ; 
printf ( " \nBearing is = $.2f Deg " , T_BEAR ) ; 



66 



y********************************************-*-*-*-*-*-*-**-*-**--*-**--*-*-*-*-*-*-*-* 



This function will calculate the range and angular error 
as in the formula in page 38. 



float tr_rang_error ( tr_bandwid, tr_freq ) 
float tr_bandwid, tr_freq; 

{ 

float x, RAN, TER1 , TER2 , TER, TSN ; 
float random(), tr_sig_noise( ) ; 
if ( TRACK ==2 ) 

{ 

RAN = random( ) ; 

TSN = tr_sig_noise( ) ; 
x = PI * tr_bandwid / tr_freq; 

TER1 = PI * tr_bandwid * pow( ( 2 . 0*TSN ) , 0 . 5 ) ; 

TER2 = sin( x*PI/ 180 ) / x; 

TER = RAN * 1.73205 * TER2 * C_SPEED / ( TER1 * 2.0); 
return( TER ) ; 

> 
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float tr_el v_error ( tr_bel ) 
float tr_bel; 

{ 

float RAN, TELE, TSN ; 

float random(), tr_sig_noise( ) ; 

if ( TRACK == 2 ) 

{ 

RAN = random( ) ; 

TSN = tr_sig_noise( ) ; 

TELE = (RAN * 0.628 * tr_bel) / (2.0 * pow( ( 2 . 0*TSN ) , 0 . 5 ) ) ; 
return ( TELE ) ; 

) 

} 

float tr_bear_error( tr_bear ) 
float tr_bear; 

{ 

float RAN, TBEA , TSN; 

float random(), tr_sig_noise( ) ; 

if ( TRACK == 2 ) 

{ 

RAN = random( ) ; 

TSN = tr_sig_noise( ) ; 

TBEA = (RAN * 0.628 * tr_bear) / (2.0 * pow((2.0 * TSN), 0.5)); 
return( TBEA ) ; 

} 
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/•M*************************************************************** 

rad_head . h 

Define the parameter to be used in the search radar and 
tracking radar programs. 



/^define 


C_SPEED 


3. 0e8 


/* 


meter/ sec 


**/ 


^define 


K_BOLTZ 


1 . 38e-23 


/* 


joules/degree 


*/ 


//define 


TRUE 


1 








//define 


FALSE 


0 








//define 


DETECTED 


1 








//define 


UNDETECTED 0 
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APPENDIX C 



USERS’ GUIDE 

This appendix provides users’ instructions for the radar detection 
model . 

DATA function is used to read the input parameters in the 
computer. The function has two beneficial features. It assigns 
the initial values to the parameters, and it allows the player to 
change the values at the beginning of the program. All values 
used in the first case of a job are continued into the second 
case, except those specifically changed. 

Table 1 shows the input parameters to consider. 

TABLE 1 

RADAR DETECTION MODEL INPUTS 
DEFINITION 

Average transmitter power (Watts) 

Receiver noise temperature (290 K) 

Pulse repetition frequency (Pulses) 

Pulse width (sec) 

Receiver recovery time (sec) 

Antenna height (feet) 

Antenna area (Square meters) 

Radar frequency (hertz) 

Radar receiver bandwidth (hertz) 



INPUT 

r_power 

r_temp 

r_prf 

r_tau 

r_taur 

r_ant 

r_area 

r_f req 

r bandwid 
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INPUT 



DEFINITION 



r_f igijre 


Receiver noise figure 


c _signal 


Critical signal to noise ratio 


r_false 


False alarm number ( 1 . e6 to 1.e12) 


r_scan 


Number of scans 


r azr 


Antenna scan rate (deg/sec) 


r_baz 


Antenna azimuth beam width (deg) 


r_bel 


Antenna elevation beam width (deg) 


r_plum 


Plumbing loss (Db) 


r_beam 


Beam shape loss (Db) 


j_power 


Jammer power (Watts /GHz) 


g_ jammer 


Jammer antenna gain (Db) 


j_range 


Range from jammer to the radar (N.M) 


tr_power 


Tracking radar transmitter power(Watts) 


tr_ant 


Antenna heigh t( feet ) 


tr area 


Antenna area( Sq . meters ) 


tr_f req 


Tracking radar frequency ( hertz ) 


tr_bandwid 


Tracking radar bandwidth (hertz) 


tr_ s ignal 


Critical signal to noise ratio 


tr_gate 


Track gate losses (Db) 


tr_beam 


Beam shape losses (Db) 


tr_bel 


Antenna elevation beam width (deg) 


tr_bear 


Antenna azimuth beam width (deg) 
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GIVEN BV THE TARGET 



INPUT 


DEFINITION 




t_rcs 


Target radar cross section (Sq. meters) 


x_t 


X-position of the 


target(Meters ) 


y_t 


Y-position of the 


target(Meters ) 


z _t 


Z-position of the 


target ( Meters ) 
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APPENDIX D 



SAMPLE INPUT VALUES 

To provide an illustration of the use of the computer model 
let us suppose that the radar system possesses the following para 
meter values 



INPUT 


VALUES 




r power 


2 . 0e5 


Watt 


r_temp 


290 


K 


r_prf 


1500 


Pulses/sec 


r_tau 


2 . 0e-6 


sec 


r_taur 


0 . 5e-6 


sec 


r ant 


30.0 


feet 


r_area 


0.8 


Square Meters 


r__f req 


9 . 0e9 


Hz 


r bandwid 


5 . 0e6 


Hz 


r_f igure 


5.0 




c_signal 


5.0 




r_false 


1 . 0e1 0 




r_scan 


3.0 


Scans 


r_azr 


120 


deg/sec 


r_baz 


1 .5 


deg 


r_bel 


10.0 


deg 


r_plum 


2.0 


Db 


r_beam 


3.0 


Db 


j__power 


10 to 0. 


1 Watts/GHz 



73 



INPUT 



VALUE 



g_ jammer 


1 . 0 


Db 


j range 


60.0 


N.M 


tr_power 


1 . 2e5 


Watts 


tr freq 


4 . 0e9 


Hz 


tr_ant 


30.0 


feet 


tr_area 


0.6 


Sq . meters 


tr_bandwid 


4 . 0e4 


Hz 


tr_signal 


5.0 




tr_gate 


1 . 0 to 1 . 6 Db 


tr_beam 


0.5 


Db 


tr_bear 


1 . 1 


Deg 


tr bel 


0.7 


Deg 



GIVEN BY THE TARGET 



t_rcs 


5.0 


Sq . meters 


x_t 


40000 


Meters 


y_t 


55000 


Meters 


z t 


5000 


Meters 
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